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Abstract 

The basic aspects of both Bohzmann-Gibbs (BG) and nonextensive statisti- 
cal mechanics can be seen through three different stages. First, the proposal 
of an entropic functional {Sbg = —kJ2iPi^^Pi for the BG formalism) with 
the appropriate constraints {J2i Pi = ^ and J2i Pi^^i = U for the BG canonical 
ensemble). Second, through optimization, the equilibrium or stationary-state 
distribution {pi = e~^^^ /Z^g with Z^g = J2j e~^^^ for BG). Third, the con- 
nection to thermodynamics (e.g., Fbg = —'^^'clZbg and Ubg = —-§p^'^Zbg)- 
Assuming temperature fluctuations, Beck and Cohen recently proposed a gen- 
eralized Boltzmann factor B{E) = d(3f{f3)e~^^. This corresponds to the 
second stage above described. In this letter we solve the corresponding first 
stage, i.e., we present an entropic functional and its associated constraints 
which lead precisely to B{E). We illustrate with all six admissible examples 
given by Beck and Cohen. 
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The foundations of statistical mechanics and thermodynamics, is a fascinating and subtle 
matter. Its far reaching consequences have attracted deep attention since more than one 
century ago (see, for instance, Einstein's remark on the Boltzmann principle [0]). The field 
remains open to new proposals focusing nonequilibrium (and metaequilibrium) stationary 
states. One of such proposals is nonextensive statistical mechanics, advanced in 1988 
(see for reviews). This formalism is based on an entropic index q (which recovers usual 
statistical mechanics for g = 1), and has been applied to a variety of systems, covering 
certain classes of both (meta)equilibrium and nonequilibrium phenomena, e.g., turbulence 
0, hadronic jets produced by electron-positron annihilation 0, cosmic rays 0, motion of 
Hydra viridissima one-dimensional and two-dimensional [1^ maps, among others. In 
addition to this, it has been advanced that it could be appropriate for handling some aspects 
of long-range interacting Hamiltonian systems (see [|lT| and references therein). 

Let us be more precise. Nonextensive statistical mechanics is based on the entropic form 
(here written for W discrete events) 

S^ = k (5:p, = l;gG7^) (1) 

^ i=i 

with 

w 

Si = Sbg = -k^pMpi . (2) 

If we focus on the canonical ensemble (system in contact with a thermostat), we must add 
the following constraint 

^'t ,' =U, {U, = Ubg), (3) 

where {Ei} is the set of eigenvalues of the Hamiltonian with given boundary conditions. 
Optimizing 5*^ we straightforwardly obtain the distribution corresponding to the equilibrium, 
metaequilibrium or stationary state, namely 



with 

w 



Z, = Y.\^-{l-qmE,-U,)Y/^'-'^^ (5) 



and 



f^'i = ^^W—o , (6) 

(5 being the Lagrange parameter. We easily verify that, for g = 1, we recover the celebrated 
BG equilibrium distribution 

P. = ^ (7) 

^BG 

with 

w 

Zbg = E ■ (8) 



Very recently, Beck and Cohen have proposed [|T^] a generalization of the BG factor. More 
precisely, assuming that the inverse temperature (3 might itself be a stochastic variable, they 
advance 

PCX) 

B{E)= / rf/?7(/?')e-^ ^ , (9) 
Jo 

where the distribution f{P') satisfies 

oo 

dp'fiP') = 1 . (10) 

-oo 

It is clear that f{f3') = 6{(3' — j3) recovers the usual BG factor. They have also shown 
that, if /(/?') is the Gamma (or x^) distribution, then the distribution associated with 
nonextensive statistical mechanics is reobtained. They have also illustrated their proposal 
with the uniform, bimodal, log-normal and F— distributions. Moreover, they define (see 
also H) 

^Bc = ((...)= j^J(3fm-))- (11) 



where we have introduced the notation qbc in order to avoid confusion with the present 
q. Clearly, if /(/?') is the Gamma distribution, then qsc = Q- Finally they argue that, 
whenever {qbc ~M << 1? the B{E) factor asymptotically becomes the nonextensive one for 
all admissible f{(3'). In other words, nonextensive statistical mechanics would correspond, 
for this particular mechanism where nonextensivity is driven by the fluctuations of /?, to the 
universal behavior whenever the fluctuations are relatively small. 

This is no doubt a very deep and interesting result, hut it does not constitute by itself 
a statistical mechanics. The reason is that the factor B{E) has been introduced through 
what can, in some sense, be considered as an ad hoc procedure. The basic element which is 
missing in order to be legitimate to speak of a statistical-mechanical formalism is to be able 
to derive the factor B{E) from an entropic functional with concrete constraints (and very 
especially the energetic constraint, which generates the concept of thermostat temperature). 
The purpose of the present paper is to exhibit such entropic form and constraint. 

Let us first write a quite generic entropic form (from now on /c = 1 for simplicity) , namely 

w 

S = Y: s{p.) {s{x) > 0; s{0) = s{l) = 0) . (12) 

i=l 

For the BG entropy Sbg we have s{x) = —xlnx, and for the nonextensive one Sq we have 
s{x) = {x — x'^)/{q — 1). The function s{x) should generically have a definite concavity 
Va; G [0, 1]. Conditions (12) imply that S* > and that certainty corresponds to S* = 0. 

Let us now address the constraint associated with the energy. We consider the following 
form: 

^hv'^^t^f =U (0 < u{x) < 1; ^.(0) = 0; «(1) = 1) . (13) 

For the BG internal energy Ubg we have u{x) = x, and for the nonextensive one Uq we 
have u{x) = x'^. The function u{x) should generically be a monotonically increasing one. 
Certainty about Ej implies U = Ej. The quantity u{pi) / J2fLiu{pj) constitutes itself a 
probability distribution (which generalizes the escort distribution defined in ||14|]). 
Let us consider now the functional 



where a and /5 are Lagrange parameters. The condition d^/dpj = imphes 

s'ip,) -a- \ , u'ip,){E, -U)=0. (15) 
Let us now heuristically assume 

u{x) = ^ + I's'^x) , (16) 

u{x) = jjix + us{x) + ^ . (17) 

But the condition s(0) = 'u(O) = imphes ^ = 0, and the conditions s(l) = and u{l) = 1 
imply that /i = 1, hence 

u{x) = X + z/s(x) (18) 

hence, 

u{x) = l + us'{x) . (19) 

Observe that if z/ = we have ^(a:;) = a; and Yl^i u{pj) = Z^JLiPj = 1- 
The replacement of Eq. (19) into Eq. (15) yields 



a + /3- 



E,-U 



^'(P.) = 1 J-lC^ . (20) 



hence 



a + /? 



E,-U \ 



l-/5z/^^^ 



(21) 



The condition Yl^i Pi = 1 enables the (analytical or numerical) elimination of the Lagrange 
parameter a. The function (21) is to be identified with B{E)/ Jl^ dE'B(E') from [|1^. In 
other words, if E{y) is the inverse function of B{E)/ dE'B{E') we have that 
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« + 



E(y)-U 



s{x) = / dy 



(22) 







E(y)-U 



which, together with Eq. (18), completely solves the problem once v is determined. Sum- 
marizing, given an admissible function B{x), we have uniquely determined the functions 
s{x) and u{x), which, replaced into Eqs. (12) and (13), concludes the formulation of the 
statistical mechanics associated with the Beck-Cohen superstatistics. An important remark 
remains to be made, namely any monotonic function of S given by Eq. (12) also is a solution 
at this stage. Which of those is to be retained for a possible connection with thermody- 
namics is a different matter, and remains an open issue at the present stage. For example, 
for nonextensive statistical mechanics, in what concerns the stationary distribution, Sq and 
the Renyi entropy — ^t^[^+^-^)s<i\ equivalent. In other words, at this level we could 
indistinctively use Sq or S^. There is however, a variety of physical arguments which are 
out of the scope of the present work but which nevertheless point, in that particular case, Sq 
as being the correct physical quantity to be used for thermodynamic and dynamic purposes. 

Let us now compare the B{E) factor obtained by Beck-Cohen formahsm with the distri- 
bution presented here in Eq. (21). In addition to the fact that the former does not include 
the normalization constant whereas the latter does, we notice that the B{E) factor has 
parameters such as /3o = instead of the parameters a, C/, v and P/Y^Y '^(Pj) (general- 
ization of (6)) appearing in Eq. (21). This problem will be handled as follows. Since our 
aim is to determine the functional forms of s{x) and u{x), it is enough to work with only 
one variable. So, we take /3o = /3/ Ef uipj) = 1 and C/ = 0. Let us now determine u. Using 
Eqs. (15) and (19), and integrating, we obtain 



It is physically reasonable to assume that u{x) monotonically increases with x, hence 
du/dx > 0, hence (1 -|- oiu)/(l — vE) > 0. We shall verify later that 1 + au > 0, hence it 
must he u < 1/E. If we note E* the lowest admissible value of E, we are allowed to consider 
1/ — 1/E*. In particular, if — oo then it must be = 0. An example where E* is 




(23) 
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finite is nonextensive statistical mechanics with g > 1. In this case, E* = 1/(1 — g), hence 
u = 1 — q. We can trivially verify that this value for u, together with s{x) = (x — x'^) /(g — 1) 
and u{x) = x"^ precisely satisfy Eq. (18). 

Summarizing, the final form of u{x) is given by 

u(.)^il^a/E-)l'^-J^, (24) 

and therefore 

, , /-^ , a + E(y) 

s(x) = / dy ; . (25) 

^ ' Jo ^1-E{y)/E* ^ ' 

In what follows we shall illustrate the above procedure by addressing all the admissible 
examples appearing in |12|. The cases associated with the Dirac-delta and the Gamma distri- 



butions for f{l3) (respectively corresponding to BG and nonextensive statistical mechanics) 
can be handled analytically. The other four cases (uniform, bimodal, log-normal and F- 
distributions for /(/?)) have been treated numerically as follows. We first choose /(/5), then 
calculate B{E) and from this calculate /q°° dE'B{E'). By inverting the axes of the variables, 
we find the inverse E{y) of Beck-Cohen superstatistics. From this we obtain E*. Two cases 
are possible. The first one corresponds to E* — oo, hence z/ = 0, hence u{x) = x and 
s{x) = ax + dyE{y). The condition s(l) = determines a, which is therefore given by 
a = — Jq dyE(y). The second case corresponds to a finite and known value of E*, which 
determines u = 1/E*. From this, we calculate Jq jz^j^Je^- From the condition u{l) = 1 
and using Eq. (24), it is easy to see that \ + av = 1 + a/ E* = 1/ i^^^^^^yE* ^ ^ ^"^^ 
is direct to determine a, which in turn enables the calculation of u{x) using once again Eq. 
(24). Finally, from Eq. (18), we calculate s(x), and the problem is solved. 

In Fig. 1 ((a) and (b)) typical examples of s{x) and u{x) are presented for all the 
cases addressed in []12[. In Fig. 2 we show the entropies associated with all these examples 



assuming W = 2. 

Let us conclude by saying that it has been possible to find expressions for the entropy 
and for the energetic constraint that lead to a generic Beck-Cohen superstatistics. Of course. 



similar considerations are valid for other constraints if we were focusing on say the grand- 
canonical ensembles. The step we have discussed is necessary for having the statistical 
mechanics generating these superstatistics through a variational principle. What remains to 
be done is the possible connection with thermodynamics. This is not a trivial task, because, 
unless we are dealing with a nonlinear power-law for u{x) (which precisely is nonextensive 
statistical mechanics), the Lagrange parameter a is not factorizable in Eq. (21), hence 
no partition function can be defined in the usual sense, i.e., a partition function which 
depends on P (and other analogous parameters), but does not depend on a. Summarizing, 



nonextensive statistical mechanics not only paradigmatically represents, as shown in [12 



the universal behavior of all Beck-Cohen superstatistics in the limit qbc — q — 1, but it is 
the only one for which an a-independent partition function can be defined. 

Last but not least, let us emphasize that the present results strengthen the idea that 
the statistical mechanical methods can be in principle used out of equilibrium as well. To 
be more specific, we can think of using them (i) in equilibrium (e.g., in the t —>■ oo limit 
of noninteracting or short-range interacting Hamiltonians, as well as in the lim^v^oo hnit^oo 
of long-range interacting many-body Hamiltonian systems; this is essentially BG statistical 
mechanics), (ii) in metaequilibrium (e.g., in the limt^oo hmTv^oo of long-range interacting 



many-body Hamiltonian systems; see, for instance, ||rT|), and (iii) for appropriate classes of 
stationary states. 
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Fig. 1 - Functions s{x) (a) and u{x) (b) for illustrative examples of the cases focused 
on in [|12[. From top to bottom: (i) F-distribution [f{(3) = 12^^]', (ii) Bimodal [/(/5) = 
0.5 6{(3 - 1/2) + 0.5 6{l3 - 3/2)]; (iii) Log-normal [/(/3) = e"^^^^^^]; (iv) Uniform 



(3V2n 

IfiP) = 1 on the interval [1/2, 3/2], and zero otherwise]; (v)Dirac delta [/(/5) = 6{P — 1)]; 
(vi) Gamma [/(/3) = (q sj^^^rii 25) ^""^'^^ ^] ; which corresponds to g = 1.8. 
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Fig. 2 - W — 2 entropies associated with the examples presented in Fig. 1. 
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